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TWO-DIMENSIONAL COMPRESSIBLE FLOW IN TURBOMACHINES WITH CONIC FLOW 

SURFACES 

Bv John D. Stanitz 


SUMMARY 

A general method of analysis is developed for two-dimensional, 
steady , compressible flow in stators or rotors of radial - and 
mixed-flow turbomachines with conic flow surfaces {surfaces of 
right circular cones generated by center line of flow passage in 
the axial-radial plane). The variables taken into account are: 
(1) tip speed of the rotor , {2) flow rate , (8) blade shape , (f) varia- 
tion in passage height with radius , {5) number of blades , and 
(6) cone angle of the flow surface. Relaxation methods are used 
to solve the. nonlinear differential equation for the stream function. 

The analysis indicates that: ( 1 ) The solution obtained for a 
given turbomachine also applies to certain other ( equivalent ) 
turbomachines with a larger or smaller number of like passages 
{same spacing of the blades on the conic flow surface , same blade- 
thickness distribution , and so forth) but with different cone 
angles; {2) for the same number of similar blades , the blade 
loading is less for mixed-flow than for radial-flow turbomachines; 
and {8) any solution obtained for an outflow turbomachine with 
shockless {smooth) entry is also the solution for an inflow turbo- 
machine with shockless entry and with the flow direction and 
blade rotation {if any) reversed. 

Two numerical examples are presented; one for compressible 
and the other for incompressible flow in a centrifugal compressor 
with thin , straight blades. The solutions were obtained in a 
region of the compressor , including the impeller tip , that was 
assumed to be unafected by the inlet configuration of the im- 
peller or by the diffuser vanes {if any) . Both examples are for 
the same impeller ( 18 ° included angle between blades on conic 
flow surface) with the same tip speed {equivalent to a tip Mach 
number of 1.5 for the compressible-flow example ), with the same 
flow rate , and with a constant flow area normal to the flow 
surface. The results of these examples are given by plots of the 
streamlines , constant velocity-ratio lines , and constant pressure- 
ratio lines. 

It is concluded from the examples that , if the fluid in high- 
speed, rotating, radial- and mixed-flow blade systems is com- 
pressible, incompressible solutions give poor quantitative results 
{exception, the slip factor) and , in som.e respects , poor qualita- 
tive results. 

INTRODUCTION 

Increased knowledge of flow conditions within radial- and 
mixed-flow compressors and turbines should indicate means 
of improving performance of these turbomachines. For 
example, boundary-layer separation, which decreases the 
efficiency of these machines, can be minimized or eliminated 


by aerodynamic design based on knowledge of the velocity 
gradients that result from various design configurations. 

For a given set of operating conditions, the flow conditions 
within radial- and mixed-flow turbomachines depend on the 
geometry of the machine (three-dimensional-flow effects) 
and on the properties of the fluid (compressibility and vis- 
cosity). Most treatments of the problem have been con- 
cerned with the two-dimensional-flow effects for incom- 
pressible, nonviscous fluids. (For example, see references 1 
to 5.) 

In the analysis reported herein, compressibility is consid- 
ered. This consideration is especially important in radial- 
and mixed-flow turbomachines because the large pressure 
ratios per stage result in density changes that greatly affect 
the fluid velocities, and so forth. The analysis is developed 
for two-dimensional, compressible, nonviscous, steady flow 
through stators or rotors of radial- and mixed-flow turbo- 
machines in which the center line of the flow passage in the 
axial-radial plane generates the surface of a right circular 
cone when rotated about the axis of the machine. The two- 
dimensional-flow pattern is considered to lie upon this surface. 

The solution of two-dimensional, compressible-flow equa- 
tions (‘an be accomplished by relaxation methods, which 
were developed by Southwell (references 6 and 7) and 
which have been applied to compressible-flow problems by 
Emmons (reference 8). It is essentially the procedure out- 
lined in reference 8 that is employed in the numerical solu- 
tion of the differential equation obtained in this analysis. 

The analysis is developed for turbomachines with arbi- 
trary blade shapes and is applied, in the numerical examples, 
to a flow region, including the impeller tip, of a centrifugal 
compressor with straight, thin blades that lie on conic radii 
(elements). A simplified analysis for straight blades lying 
on conic radii is developed that checks the results of the 
relaxation solution within the impeller except for the flow 
region near the impeller tip. 

This analysis was developed at the NACA Cleveland 
laboratory in 1947. 

ANALYSIS 

PRELIMINARY CONSIDERATIONS 

This analysis develops a general method whereby the 
streamlines, velocity distribution, and pressure distribution 
can be determined for steady, two-dimensional compressible 
flow in stators or rotors of radial- and mixed-flow turbo- 
macliines with arbitrary blade shapes and varying passage 
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heights. The radial component of the flow may be in the 
direction of increasing radius (outflow turbomachine) or in 
the direction of decreasing radius (inflow turbomachine). 

Conic flow surface. — The analysis is limited to turboma- 
chines in which the center line of the passage in the axial- 
radial plane generates the surface of a right circular cone, 
with the cone angle a (fig. 1), when rotated about the axis 
of the machine. (All symbols are defined in appendix A.) 
The two-dimensional-flow pattern is considered to lie upon 
this cone surface (hereinafter referred to as “conic flow sur- 
face”)- For mixed-flow turbomachines the cone angle a is 
less than 180° but greater than 0°. For the special case in 
which a is 180°, the conic flow surface becomes flat and is 
normal to the axis of the machine. Such turbomachines 
(a=180°) are designated radial-flow machines. For the 
special case in which a is 0°, the conic flow surface becomes 
cylindrical and is concentric with the axis of the machine. 
Such turbomachines (a=0°) are designated axial-flow 
machines. Axial-flow machines are not considered in this 
analysis for reasons that are subsequently discussed. 

Coordinate system. — The developed view of a conic flow 
surface is shown in figure 2. The dimensionless, conic 
coordinates R and 0 of this conic flow surface are relative to 



Figure 1.— Fluid particle on coordinate system relative to blades. Blades may be stationary 
(stator blades) or rotating (rotor blades). Center line of flow passage generates surface of 
right circular cone with cone angle a. 


Direction of rotation 
for rotor bfoefes 




Figure 2.— -Fluid particle on developed view of conic flow surface. This surface may be 
stationary (stator blades) or rotating (rotor blades) . R and 0 , dimensionless coordinates 
relative to blades; H, passage-height ratio normal to conic flow surface; l/and V, tangential 
and radial components, respectively, of relative velocity ratio Q. 

the blades. The blades, and therefore the coordinate sys- 
tem, may be stationary (stator blades) or rotating (rotor 
blades). The conic-radius ratio R is defined as 

R=- (1) 

r T K J 

where r is the conic radius (distance along conic element from 
apex of cone) and where the subscript T refers to the blade tip 
(either the nose or the tail of the blade, whichever has the 
larger conic radius). The passage-height ratio H (figs. 1 and 
2) is normal to the conic flow surface and is a continuous 
function of the conic-radius ratio R 

ff=£=/(S) (2) 

where h is the passage height at any conic-radius ratio R. 

Velocity ratios. — A fluid particle on a developed conic 
flow surface has a relative tangential-velocity ratio U (fig. 2) 
and a radial-velocity ratio V (fig. 2) that are related to the 
relative velocity ratio Q by 

Q=± = {U*+V>)» (3a) 

o 

where 

t/= c - (3b) 

V=j (3c) 

V o 

where 

c local speed of sound 

q velocity of fluid relative to blades 

u tangential component of q (positive in direction of 

increasing 6) 

v radial (along conic element) component of q (positive 
in direction of increasing conic-radius ratio) 
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Subscript : 

o absolute stagnation condition in region of uniform 
flow upstream of blades 

The relative velocity ratio Q is defined on a coordinate sys- 
tem that is relative to the blades, therefore the velocity g 
(and u) is absolute for a stationary coordinate system (stator 
blades) and relative for a rotating coordinate system (rotor 
blades) . 

Assumptions and limitations. — This analysis assumes that 
the flow varies only along the conic flow surface, that is, 
that flow conditions are uniform across the passage normal 
to the conic flow surface. In order to satisfy this assump- 
tion it is necessary that: (1) The gradient of h with respect 
to r be small; and (2) the cone angle a (fig. 1) be sufficiently 
large. The allowable variation in a from 180° will depend 
on the relative magnitudes of h and r and on the desired 
accuracy. For small values of a the flow must be assumed to 
exist in concentric annuli, each with negligible passage 
height. For the hypothetical limiting case in which the 
ratio hjr approaches zero everywhere along the conic flow 
surface, the analysis is accurate for all values of a. 

The analysis assumes that steady flow exists relative to 
the blades. The relative motion between stator and rotor 
blades introduces pulsations that make the flow unsteady. 
These pulsations rapidly diminish, however, as the stators 
and the rotors are moved apart so that the relative flow may 
be treated as steady (between boundaries, which are far 
enough upstream and downstream of the blades to obtain 
uniform flow) provided the stators and the rotors are not 
too close together. 


dT = 0 = \ {TUSd'p-\- U') Rd6]^dR — \VdR\gd6 


where the blade-tip Mach number M r is defined by 




co r T sm g 
c 0 


( 5 ) 


where co is the angular velocity of the rotor and where 
(RM t + U) is the tangential component of the absolute 
velocity ratio. (For stator blades M T is zero.) After simpli- 
fication, 

-2M t =^+U b -^- ( 6 ) 


Substitution of the stream function i p as defined by equa- 
tions (4a) and (4b) gives 


2M t H £ ^ - tn (loge H) * - 


fc ( log -£).-»K log ‘£l (7) 


where the double coordinate subscripts ( RR and 66, in this 
case) refer to second partial derivatives with respect to the 
coordinates. 

General energy equation. — The general energy equation is 
used to determine the density ratio p/p 0 in the differential 
equation (7). When expressed in terms of the velocity 
ratios defined by equations (3b) and (3c), the general energy 
equation becomes 


DIFFERENTIAL EQUATIONS FOR FLOW IN Rd- PLANE 

The differential equations for steady, two-dimensional, 
compressible flow are developed from the continuity equa- 
tion, the equation for absolute irrotational motion, and the 
general energy equation. 

Continuity and stream function. — From steady-flow con- 
tinuity considerations for the fluid particle in figure 2 

(- p - VHr) +(-- UH ) =0 (4) 

\Po / R \Po /Q 

where p is the weight density of the fluid, and where the 
coordinate subscripts ( R and 6, in this case) refer to partial 
derivatives with respect to the coordinates. 

A dimensionless stream function satisfies the continuity 
equation (4) if defined as 

^e= £ VHR (4a) 

Po 

and 

(4b) 

Po 

Absolute irrotational motion. — In the absence of viscosity, 
shock, nonuniform heat addition, and so forth, the absolute 
motion of a fluid particle is irrotational. The dimension- 
less absolute circulation dr about the particle in figure 2 
is therefore zero, and 


Jc„T+ c / g [(RM r + UY+ V*)=Jc v T 0 + Mt ^ (X -\ v ) (8) 


where 

J mechanical equivalent of heat 
c v specific heat at constant pressure 
T static (stream) temperature 
g gravitational acceleration 
Subscript: 

U upstream boundary (boundary in region of uniform flow 
upstream of blades) 

and where the “whirl” ratio X is defined by 

\=R(RM r -{- U) (9) 


which is the whirl or absolute moment of momentum ^radius 
times absolute tangential velocity, r sin-^X^r sin 
divided by a constant (r T sin The last term in equa- 

tion (8) is the work done on the fluid and is equal to ^ tC ° 


times the change in whirl ratio. The total work done on 
the fluid is given by the last term in equation (8) with X equal 
to \ D (where subscript D refers to downstream boundary, 
the boundary in the region of uniform flow downstream of the 
blades). This total work is positive for compressors and 



4 


REPORT 935 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


negative for turbines. The whirl ratios \ v and \ D are con- 
stant in the uniform flow regions upstream and downstream 
of the blades (constant absolute moment of momentum). 
‘Rearrangement of equation (8) with 

c 0 2 =(y—l)Jgc p T 0 

where y is the ratio of specific heats, results in 


given blade shape, determines £ and 77 as functions of R 
and 6 


y=y(R,o)) 


(13a) 


Equation (7), in terms of the transformed coordinates £ 
and 7] given by equation (13a), becomes (appendix B) 


£-l+V [(KM T y-Q*-2M T \ v ] 


(10) 


2 M t H p 

q_i Po 


'Pan+'l'm—'I'i 


( lo <;)r*'( log *£>'- 


from which 


P 

Po 


<sr-i 


l 2 [( RMt ) 2 Q 2 2ibf7’Xc/]| (11) 


(fov i pu l ) ^ loge ^/)^._(io ge H) v u r ] 


and equation (12) becomes 


(14) 


Also, from equations (3a), (4a), and (4b) 

Q p» = [(^) + G&) ] (12) 

Equations (11) and (12) together with the general differential 
equation (7) provide three equations with three unknowns: 
ip, Q, and p/p 0 . The solution of these equations determines 
the steady, two-dimensional flow of compressible fluid 
through turbomachines with arbitrary blade shape, with 
arbitrary variation in the passage-height ratio, and with 
constant cone angle. 


METHOD OF SOLUTION 

Equation (7), which is nonlinear, can be solved (together 
with equations (11) and (12)) by relaxation methods. 

Relaxation methods. — Values of i p are estimated at each 
point of a grid system placed within the boundaries of the 
problem, and the residuals R, which result from the estimated 
values of \p, are computed for each grid point by expressing 
the differential equation for \p in finite-difference form with 
the sum of all terms equal to R instead of zero. The solu- 
tion is then obtained by systematically varying (relaxing) 
the values of \p at the grid points inside the boundaries 
until the values of R approach zero. 

Transformation of coordinates. — For the numerical solu- 
tion of this problem by relaxation methods, it is convenient 
(but not necessary) to select a new set of coordinates (refer- 
ence 8) so that blades of arbitrary shape in the physical plane 
( R,6 coordinates) become thin, straight, and parallel in the 
transformed plane (£,77 coordinates). Thus a grid of equally 
spaced points can be placed between the blades. This trans- 
formation of coordinates is represented by the general 
analytic function 

£(i?, 0) +irj(R, 6) =j[R exp (i0)] (13) 

where the Cartesian coordinates £ and 77 in the ^-plane cor- 
respond to velocity potential lines (£=constant) and stream 
lines (7?= constant) in the 7?0-plane for incompressible flow 
past the blades, which, for purposes of the transformation, 
are considered to be stationary (co=0) and to have a con- 
stant height (H= 1). Equation (13), in specific form for a 


Qj = %W + +rh* (15) 


where H is now a function of £ and 77 (given by equations (2) 
and (13a)) and where the coefficients u i} v u and q t are deriv- 
atives of equation (13a) defined by 


Ui — ^ R ~~R 

(16a) 

II 

1 

(16b) 

g_i={u l ^+v i 2 y- < ‘ 

(16c) 


where the subscript i indicates that the coefficients corre- 
spond to incompressible velocities in the 7?0-plane. 

For certain simple blade shapes, equation (13) is a simple 
analytic expression that determines £(#,0) and 77 (Rfi) (equa- 
tions (13a)) directly. For arbitrary blade shapes, however, a 
specific expression for equation (13) is not readily available 
and it is easier to obtain £(/?,0) and r)(R, 9) by relaxation 
solutions of the Laplace equations for £ and 77 in the Rd- 
plane (appendix C). 

Finite-difference equations. — In order to solve the system 
of equations (equations (11), (14), and (15)) by relaxation 
methods, equations (14) and (15) must first be changed to 
finite-difference form. This change is accomplished with 
the aid of the following equations (reference 7, p. 19), which 
are based on first-order differences: (Note that higher-order 
differences could be used, which would result in more com- 
plex finite-difference equations but which would enable 
larger grid spacing, and therefore fewer grid points, for the 
same degree of approximation.) 


r <~h w-™ 

(K + F.-2F) 


( 17 ) 
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Figure 3.— Sample grid showing grid spacing b and numerical subscript convention for 
adjacent grid points. 


Q f= 26 B Mt-W (19) 

After the values of \f/ have been estimated at the grid 
points inside the boundaries, the system of equations (11), 
(18), and (19) provides three equations with three unknowns 

—> Q, and R at each grid point. Equations (11) and (19) 

Po 

determine the values of the density ratios in equation (18), 
which, is then solved for the residual R . 

BOUNDARY CONSIDERATIONS 

The values of \p at the grid points inside the boundaries 
depend upon the values of i {/ at the grid points along the 
boundaries. These boundary values of \f/ are determined 
by the design characteristics and the operating conditions 
of the turbomachine. 

Location of boundaries. — The boundaries of the flow field 
in the i?<9-plane (fig. 4(a)) are the blade surfaces and the 
upstream and downstream boundaries at constant values of 
R , which are any distance far enough from the blades to 
insure uniform flow conditions at these boundaries. The 
upstream and downstream boundaries enclose all the blades; 


where 

F any twice-differentiable function of two variables (£ and 
r) f in this case) 
b grid spacing 

Subscripts: 

1,2, 3,4 grid points adjacent to point being considered 
(F with no subscript) 

A sample grid is shown in figure 3. The grid spacing b is 
arbitrary. However, the smaller the value of b , that is, the 
larger the number of grid points, the greater is the accuracy 
of the approximate, finite-difference equations (17). 

With the aid of equations (17), equation (14) becomes 

+ + — 4^— 1 . -^-Aogg Pl — loge P A — 

4 : \ Po Po/ 

* 4^/ Wl — — W' 2 

[ (loge Hi— log e Hi) v t — (loge H 2 — log e H^Ui] — 

2 M T Hb 2 p p 

~^Vo~ R (18) 

where the residual R has a nonzero value when the values of 
\p do not satisfy the differential equation (14) from which 
equation (18) was obtained. 

Equation (15) in finite-difference form becomes 


Downstream boundary for outflow 
turbomachine , or upstream boun- 
dary for inflow turbomachine, 
or ft y j 



(a) Physical TVfl-piane. 
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however, from symmetry considerations, the flow conditions 
along lines of constant R are cyclic with a period equal to 
the blade spacing so that the solution need be obtained only 
in a region that encloses the equivalent flow between two 
blades. This region is bounded by any two blades, the 
upstream and downstream boundaries, and by quasi bound- 
aries between the two blades and the upstream and down- 
stream boundaries (fig. 4(a)). These quasi boundaries may 
have any reasonable, continuous shape but must have the 
same angular spacing Ad (fig. 4(a)) as the blades, where 

Ae== ~E sin f ( 20 ) 

where B is the number of blades (or passages). It is 
convenient to select as the shape of these quasi boundaries 
the incompressible, stagnation streamlines (constant 7?) that 
are determined in appendix C with stagnation points at the 
nose and the tail of the blade. 

In the ^-plane (fig. 4(b)), the blade surfaces become lines 
of constant 77 and the quasi boundaries become extensions of 
these same lines. The upstream and downstream boundaries 
in the ^77-plane are straight (appendix C) but, in general, 
are not at right angles to the lines of constant 77. The two 
blade surfaces are generally different lengths in this plane 
(fig. 4(b)). 

yfz along blade boundaries. — The boundary values of 1 p 
along the blade surfaces are constant and can be determined 
from the following considerations: 

The differential flow rate between adjacent streamlines 
is shown in figure 5 and is given by 

dw— pvh T r T HRd6— puh T r T HdR 


where w is the flow rate between streamlines. From equa- 
tions (3b), (3c), (4a), and (4b), 


or 


dw — poC<fa T r T (\pedd-\- yp R dR) 

dw—p 0 c 0 h T r T d4/ (21) 


If w and 4 / are assigned values of zero along the positive 
blade surface (the blade surface in the direction of increasing 
6), equation (21) can be integrated across the passage to the 
negative blade surface (the blade surface in the direction of 
decreasing 6) to give 

W 

~J3 = Po c oh'T r T4 / n ( 22 ) 


where W is the total flow rate through the turbomachine 
and the subscript n refers to the negative blade surface. 
Equation (22) can be simplified by the following consider- 
ations : 

The flow area a T of the annulus at the tip of the blades is 
given by 

a T =2ir sin ^ VjJbT 


from which 
becomes 


equation (22) combined with equation (20) 


where the flow coefficient 4> is defined by 


( 24 ) 


Equation (23) determines the boundary value of \p on the 
negative blade surface as a function of the operating param- 
eter <t> and the design parameter AO. 

Equation (23) was developed for through flow in the 
direction of increasing radius ratio (outflow turbomachines) , 
that is, for positive values of the radial-velocity ratio V. 
For through flow in the direction of decreasing radius ratio 
(inflow turbomachines), that is, negative values of V , the 
magnitude of 4'n is given by equation (23) but the sign is 
changed from positive to negative. 

\f/ along quasi boundaries extending from positive blade 
surface. — Because the quasi boundaries in the ^77-plane en- 
close the equivalent flow between two blades (see section 
Location of boundaries), the values of ^ at points along the 
quasi boundaries extending from the negative blade surface 
(fig. 4(b)) are 4/ n greater (outflow machine), or 4'n less 
(inflow machine) than the values of 4/ at corresponding grid 
points (corresponding to the same value of R) along the 
quasi boundaries extending from the positive blade surface. 
Therefore, the values of 4' along the quasi boundaries extend- 
ing from the negative blade surface are not recorded or 
relaxed. 

Estimated values of the stream function 4 / along the quasi 
boundaries extending from the positive blade surface in the 
£r;-plane (fig. 4(b)) can be obtained by assuming, as a first 



Figure 5. — Fluid particle between adjacent streamlines. Radial component of flow rate, 
pvh T r T HRd $ ; tangential component of flow rate, — puh T r T HdR . 


4'n=<j> AO 


(23) 


i 
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approximation, that the flow conditions upstream and down- 
stream of the blades in the i?0-plane (fig. 4(a)) are uniform, 
that is, the flow conditions are a function of only R. The 
variation in p with £ along the quasi boundaries extending 
from the positive blade surface can then be determined 
(appendix D) from continuity and from the whirl ratios 
and Xz>, which for uniform flow remain constant upstream 
and downstream of the blades (conservation of moment of 
momentum). The whirl ratio \ v is specified and the whirl 
ratio Ad is determined, for a given blade shape and operating 
condition, by the Joukowski condition, which requires that 
the rear stagnation point occur at the blade tail (or, in the 
case of infinitely thin blades or blades with cusped tails, the 
flow must be tangent to the blade surfaces at the tail). The 
value of Ad can be estimated from considerations given in 
appendix D. 

The values of i p along the quasi boundaries extending from 
the positive blade surface (which values are obtained from 
the preceding variation in p with £) are estimated values 
and must therefore be relaxed. 

p along upstream and downstream boundaries — The 
value of p at any point along the upstream or downstream 
boundary (fig. 4(b)) is determined by the integrated varia- 
tion in i p along the quasi boundary from the fixed (zero) 
value of p on the positive blade surface and by the integrated 
variation in p along the upstream or downstream boundary 
to the point in question. The variation in p along the quasi 
boundaries was estimated in the previous section, and the 
variation in p along the upstream and downstream bound- 
aries is constant (uniform flow conditions assumed at these 
boundaries) and is of such magnitude that the change in p 
from one quasi boundary to the next is equal to p n . 

The values of p along the upstream and downstream 
boundaries (i p v and p D > respectively) arc considered fixed 
during a relaxation solution. But these values of p, for the 
initial relaxation solution, are dependent upon the estimated 
variation in p along the quasi boundaries extending from the 
positive blade surface. In general, therefore, these values of 
i p v and p D do not result in a solution that exactly satisfies the 
prescribed whirl ratio Ac upstream of the blades and the 
Joukowski condition (which, together with the blade shape 
and the operating conditions, determines Ad) downstream of 
the blades. It is therefore usually necessary, after the initial 
relaxation solution, to adjust (by methods developed in 
appendixes E and F) the values of p v and p D (keeping, 
however, the same uniform variation in p along these 
boundaries). The relaxation solution is then repeated using 
these new values of p v and p D that satisfy Ac and the Jou- 
kowski condition. 


ADDITIONAL CONSIDERATIONS 

Equivalent turbomachines with different cone angles. — 
The flow field for the flow that passes between any two 
blades is the same for all blade passages in a given turbo- 
machine. Therefore, the solution obtained for the flow 
field in a given turbomachine also applies to certain other 
(equivalent) turbomachines with a larger or smaller number 
of like passages having the same angular spacing of the 
blades A 0, blade- thickness distribution, and so forth, but 


with different cone angles a. The cone angles for the equiv- 
alent turbomachines are determined by the number of 
passages B in the machine and are given by equation (20) as 


a= 2 sin 1 


BAO 

2tt 


Also, from equation (20), 


ry 27T . a. 
B= A6 Sm 2 


so that a radial-flow turbomachine (a=180°) has more 
blades than an equivalent mixed-flow turbomachine 
(a<180°), which has the same blade loading, and so forth. 
Furthermore, if the number of blades in the equivalent 
mixed-flow turbomachine is increased to equal the number 
of blades in the radial-flow turbomachine, the blade loading 
in the mixed-flow machine is decreased, so that, in general, 
for the same number of similar blades, the blade loading is 
less for mixed-flow than for radial-flow turbomachines. 

Equivalent outflow and inflow turbo machines. — Any solu- 
tion obtained for an outflow turbomachine with shockless 
(smooth) entry is also a solution for an inflow turbomachine 
with shockless entry and with the flow direction and blade 
rotation (if any) reversed. The shockless entry for the out- 
flow machine corresponds to the Joukowski condition for the 
inflow machine and, vice versa. 

Axial-flow turbomachines. — For axial- flow turbomachines, 
the cone angle a becomes zero and the flow field is assumed 
to lie on a cylindrical surface about the axis of the machine. 
For a cylindrical surface, the conic radius r is infinite and 
therefore the angle d is zero. As a result, the cylindrical 
flow surface degenerates into a single point (1,0) on the 
developed conic flow surface (R,0) in figure 2, so that no 
solution can be obtained for axial -flow turbomachines on the 
developed conic flow surface for which this analysis was 
developed. 

NUMERICAL PROCEDURE 

A detailed outline of the numerical procedures for the 
relaxation solution of compressible-flow problems is given 
in reference 8. The emphasis is placed herein on those 
features of the solution that are peculiar to the flow in 
turbomachines with conic flow surfaces. 

The complete relaxation solution is conveniently divided 
into two sections. In the first section, the initial relaxation 
solution is obtained using approximate values of p D and p v 
that are estimated to satisfy the Joukowski condition and the 
prescribed whirl ratio A c* Id the second section, the ap- 
proximate values of Pd and pu are adjusted to satisfy the 
Joukowski condition and the prescribed Ac, and the final 
relaxation solution is obtained. A brief outline of the nu- 
merical procedure for the initial relaxation solution follows. 

I— INITIAL RELAXATION SOLUTION 

Design characteristics and operating conditions. — In 

order to solve the system of equations (11), (14), and (15) 
for the stream function p , it is necessary that the following 
design characteristics and operating conditions be specified: 


863977—50 2 
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Design characteristics: 

(1) Passage-height ratio H , which is a function of the 
conic-radius ratio R 


H=m 


( 2 ) 


(2) Cone angle a, which is constant (fig. 1) 

(3) Arbitrary blade shape, which determines (appendix C) 


y=y(R, 0 ); 


(13a) 


from which the coefficients u u v i} and q t in equations 
(14) and (15) are obtained by equations (16a), (16b), 
and (16c), respectively. 

(4) Number of blades B, which together with the cone 

angle a, determines the angular blade spacing Ad 

Ad== ~B sm 2 ( 2 °) 

Operating conditions: 

(5) Whirl ratio upstream of blades X^, where X is defined 

by 

X — R {RJ\ / lj<-\- U') (9) 


The value of \ v results from the configuration of the 
turbomachine ahead of the blades (design character- 
istic) and from the flow rate through the machine 
(operating condition) . 

(6) Tip Mach number M T , which is defined as 


M t = 


cor T sin ^ 


c 


o 


For stator blades M r is zero. 

(7) Flow coefficient <£, which is defined as 


<i P — 


W 

po^T^o 


( 5 ) 


(24) 


This coefficient is proportional to the standard equiv- 
alent flow-rate parameter W-y/d/d (reference 9) 
where 

d ratio of upstream absolute stagnation temperature 
to standard sea-level temperature 


8 ratio of upstream absolute stagnation pressure to 
standard sea-level pressure 

(8) Ratio of specific heats 7, which for a given problem is 
considered constant 

Boundary values of \p . — The locations of the boundaries in 
the f?7-plane are discussed under Location of boundaries in 
the section ANALYSIS. The boundary values of are 
determined by the design characteristics and operating con- 
ditions outlined in the previous section and by the Joukowski 
condition. The various boundary values of ip are shown on 
the relaxation grid for an outflow turbomachine in the 
fry-plane in figure 6. The manner in which these boundary 
values are obtained is summarized as follows : 

(1) The value of the stream function along the positive 
blade surface in figure 6 is arbitrarily set equal to zero. 
(See section \p along blade boundaries.) 

(2) The value of the stream function along the negative 
blade surface in figure 6 is given by 

(23) 

The stream function >p n is positive for outflow turboma- 
chines and negative for inflow machines. (See section \p 

along blade boundaries.) 

(3) The values of the stream function along the quasi 
boundaries extending upstream and downstream of the posi- 
tive blade surface depend on the specified whirl ratio 
upstream of the blades and, for a given blade shape and 
operating conditions, on the Joukowski condition downstream 
of the blades. The method for estimating \p along these quasi 
boundaries is given in appendix D. Values of \p are not 
recorded or relaxed along the quasi boundaries extending 
from the negative blade surface for reasons given in the 
section \p along quasi boundaries extending from positive 
blade surface. 

(4) The values of the stream function along the upstream 
and downstream boundaries vary uniformly (steady-flow 
condition) in the direction of increasing 77 at the rate of \p n 
per unit of r\. (This rate is positive for outflow turbo- 
machines and negative for inflow machines.) The magni- 
tude of \p is fixed at the intersection of the quasi boundaries 
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(extending from the positive blade surface) with the upstream 
or downstream boundary by the variation in yp along these 
quasi boundaries. 

Grid layout. — A system of equally spaced grid points 
placed on the flow field in the £17- plane is shown in figure 6. 
The grid spacing b is selected so that an integral number of 
spacings occur between the blades. The smaller the grid 
spacing, the greater is the accuracy of the finite-difference 
equations (18) and (19), but also, the greater is the number 
of grid points and therefore the labor involved in obtaining 
the solution. 

For reasons given in the section \ p along quasi boundaries 
extending from positive blade surface, grid points are not 
located on the quasi boundaries extending from the negative 
blade surface. 

It is convenient (appendix E) to locate the grid system so 
that a grid point lies at the blade tail on the positive blade 
surface. In this case, the grid point on the blade nose of the 
positive blade surface is generally not equally spaced from 
the interior grid points. Also, because the upstream and 
downstream boundaries are not normal, in general, to the 
£-axis (fig. 6), the grid points on these boundaries are usually 
not equally spaced from the interior grid points. In order to 
account for these unequal spacings of the boundary grid 
points, the finite-difference equations (18) and (19) must be 
modified at the adjacent interior grid points (to be subse- 
quently discussed). 

After estimating (or assuming) the values of yp at the 
interior grid points, the problem resolves itself into two 
parts: (1) calculation of the residuals R at the interior grid 
points and along the quasi boundaries extending from the 
positive blade surface (which residuals result from the 
estimated values of 1 p at these points); and (2) relaxation 
(elimination) of these residuals by suitable adjustments in 
the values of \p at these grid points. 

Residuals. — The residuals at equally spaced interior points 
of the grid system are computed from equation (18). The 
density ratios in equation (18) are determined by equation (11) 
with the aid of equation (19). At the interior grid points 
adjacent to the blade surfaces, the solution of equation (18) 
requires the density ratios at the grid points on the 
blade surfaces; these ratios can be determined by extra- 
polating the values of p/p 0 obtained at interior grid points or 
by equations (11) and (19) at the boundary grid points 
using extrapolated values of 1 p beyond the boundaries. 

In general, the grid points along the upstream and down- 
stream boundaries arc unequally spaced from the adjacent 
interior grid points (fig. 6) so that at these interior grid 
points the finite-difference equations (18) and (19) must be 
modified to account for the unequal grid spacing (reference 7, 
pp. 73-74). This unequal grid spacing also exists, in 
general, between the nose of the positive blade surface and 
the adjacent grid point on the quasi boundary (fig. 6). 

In order to compute the residuals at grid points along the 
quasi boundaries extending from the positive blade surface 
(fig. 6), equations (18) and (19) require values of that lie 
outside the flow field enclosed by the quasi boundaries 
(fig. 7). From symmetry considerations and because the 
quasi boundaries enclose the equivalent flow between two 


blades (see section location of boundaries), the values of 
yp 4 equal yp n less (or greater, for inflow turbomachines) than 
the values of yp at corresponding positions (same value of 
R , that is, same increment of £ from the tail, or nose, of the 
blades) along the row of interior grid points adjacent to the 
quasi boundary extending from the negative blade surface. 
That is, . 

yp A —yp A —'Pn (25) 

where yp A , corresponding to 1 p± for a given grid point along 
the quasi boundary extending from the positive blade sur- 



Figure 7.— Relaxation grid in £»?-plane showing values of ^ used in equation (27). 


face, is shown in figure 7. In general, \p A does not lie on a 
grid point and the value of \p A is therefore obtained by linear 
interpolation between the adjacent grid points \p B and \p c 
(fig. 7). Therefore, 

y (26) 

where b' is defined in figure 7. From equations (25) and (26), 

'p4= z 'pB—'Pn+^ (jPc-Pb) (27) 

Equation (27) determines the value of yp* required by equa- 
tions (18) and (19) in order to compute the residuals along 
the quasi boundaries extending from the positive blade 
surface. 

If all the estimated values of \p are correct, the value of R 
is zero at all grid points. If, however, the estimated values 
of \p are incorrect, the values of R are finite and may be 
positive or negative. 

Relaxation. — After the residuals are computed, it remains 
to relax (that is, reduce) these residuals by suitable changes 
in the values of \p. In order to determine the magnitude 
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of the required changes in xp at the equally spaced interior 
grid points, all terms of the finite-difference equation (18) 
are assumed to remain constant except the terms xp l -\-xp 2 + ^ 3 + 
xpi—Axp. A change in the value of 1 p therefore causes a 
four-fold change of opposite sign in the value of R, and this 
change in xp also causes an equal change in the values of R 
at each of the adjacent grid points (because relative to these 
points the change in xp amounts to a change in xp h \p 2 , fa, or 
\p 4 ). At grid points that are unequally spaced from adjacent 
points (for example, at the grid points adjacent to the up- 
stream and downstream boundaries, fig. 6) a change in 
xp changes R an amount that depends on the coefficient of xp 
in the finite-difference form of equation (14) developed for 
unequal spacing. (See previous section.) Also, the result- 
ing change in R at adjacent grid points depends on the coef- 
ficients for the terms fa, fa , fa, and xp 4 in this finite-difference 
equation. In particular, it should be noted from equation (27) 
that changes in xp B and xp c have a weighted effect upon 
the residuals at the corresponding adjacent grid points along 
the quasi boundaries extending from the positive blade 
surface, and vice versa. 

These changes in ^ and R are recorded on the grid sheet as 
the work progresses. By continually reducing (relaxing) 
the larger residuals (any desired amount), the values of all 
residuals gradually approach zero. When this condition is 
reached, the residuals are recomputed using the complete 
finite-difference equation and taking into account the new 
values of the density ratio. After the new values of R have 
been computed, the relaxation procedure is repeated and this 
cycle is continued as often as necessary to achieve the de- 
sired accuracy. 

II— FINAL SOLUTION 

The whirl ratio \ n upstream of the blades and the 
Joukowski condition downstream of the blades are governed 
by the values of the stream function specified along the up- 
stream and downstream boundaries (xp v and xp D ). In 
section I, these values of xp v and xp D were determined from 
the estimated variation in xp along the quasi boundaries ex- 
tending from the positive blade surface (see section 
Boundary values of fa and, in general, do not result in a 
solution that exactly satisfies the prescribed value of and 
the Joukowski condition. In section II, xpu and xp D are there- 
fore adjusted to satisfy these conditions and the relaxation 
solution is repeated to obtain the final distribution of 1 p in 
the flow field. 

Joukowski condition. — If the Joukowski condition is satis- 
fied, the rear stagnation point occurs at the tail of the blade, 
or, in case of infinitely thin blades or blades with cusped 
tails, the flow is tangent to the blade surfaces at the tail. 
In either case, from appendix E, 

(E2) 

where fa is the value of xp after the Joukowski condition is 
satisfied and where the subscripts a, b, c, and cl refer to the 
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Figure 8. — Relaxation grid in fij-plane showing values of ^ used in equation (E2) to check 
Joukowski condition. 

grid points along the quasi boundary extending from the 
positive blade surface on the £q-plane (shown in fig. 8). If 
equation (E2) is not satisfied by the values of xp a , fa, \ p c , and 
\p d resulting from the initial relaxation solution (section I), 
the values of xp D along the downstream boundary are adjusted 
by methods given in appendix E. As a result of adjusting 
xpD, the values of xp at all other grid points in the flow field are 
changed by amounts that are estimated by methods devel- 
oped in appendix E. 

Upstream whirl ratio X^,- — If the upstream whirl ratio is 
satisfied, the whirl ratio at any point in the region of uniform 
flow upstream of the blades is equal to the prescribed value 
Xc/, and xp$ at that point is given by equation (F2) developed 
in appendix F 

where xp is the value of fa if the specified value of X^ is 
obtained. In general, equation (F2) is evaluated at 
the upstream boundary where, because conditions are 
uniform, xp% is constant. If x p$ obtained from the initial 
relaxation solution (section I) is not equal to the value xpt* 
given by equation (F2), the values of xp v along the upstream 
boundary are adjusted by methods given in appendix F. 
As a result of adjusting xp Uy the values of xp at all other 
grid points in the flow field are changed by amounts that are 
estimated by methods developed in appendix F. 

It should be noted that the corrections for \ v affect the 
Joukowski condition, and vice versa. For low-solidity 
blades these interrelations should be considered, but for 
liigh-solidity blades the effect of changes in xp D on and the 
effect of changes in xp v on the Joukowski condition arc gen- 
erally small and can be neglected. 

After the values of xp D and xp v have been adjusted and the 
resulting changes in xp at the grid points in the flow field have 
been estimated, the relaxation methods are repeated to 
eliminate the small residuals that result from the new values 
of xp at the grid points in the flow field. After the correct dis- 
tribution of xp has been determined, the pressure and velocity- 
ratio distributions can be obtained from the density 
ratio and equations (4a) and (4b). If more detailed infor- 
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mation of flow conditions in certain regions of the flow field 
is desired, the grid spacing b can be reduced and the relaxa- 
tion methods repeated in these regions. 

Accuracy. — No quantitative evaluation of the accuracy of 
relaxation solutions is available (reference 10, p. 176). 
Because the computed velocities depend on differences in 
the values of 0 at adjacent grid points, that is, the small 
difference of large numbers, however, it is important to know; 
the Values of 0 with sufficient accuracy to assure the desired 
accuracy for the velocity calculations. In the numerical 
examples of the present report, the values of 0 were computed 
to the nearest 0.00001 compared with the maximum value 
of 0, at the negative blade surface, of 0.15700. 

NUMERICAL EXAMPLES 

Two numerical examples are presented; one for compressi- 
ble and the other for incompressible flow through the im- 
peller of a centrifugal compressor. Both examples are for 
the same impeller geometry with the same tip speed and 
weight flow. 

Flow field. — A diagram of the impeller and vaneless por- 
tion of the diffuser is shown in figure 9. The cone angle a , 
shown in figure 9, is 180° (radial-flow compressor), but the 
solution applies to certain other cone angles less than 180° 
(mixed-flow compressors) given by equation (20) for an 
integral number of similar passages B with the same included 
angle Ad between blades on the conic flow surface. (See sec- 
tion Equivalent turbomachines with different cone angles.) 
The solutions are obtained in a flow field (fig. 9) that is 
considered to be unaffected by the inlet configuration of 
the impeller and by the diffuser vanes (if any) ; that is, the 
impeller inlet and the diffuser vanes must be far enough 
removed not to affect the flow appreciably in the flow field 
investigated. In this flow field, the impeller blades are 
thin and straight and the passage-height ratio II varies in 
such a manner that the flow area normal to the conic flow 
surface remains constant. 

The values of the stream function along the boundary 
between blades (R= 0.6752 in fig. 9) are determined from a 
simplified analysis (appendix G), which assumes that for 
straight thin blades the component of the relative flow 
normal to the blades is zero. This assumption is satisfactory 
(appendix G) at radius ratios within the impeller sufficiently 
far from the tip (at radius ratios less than 0.80 for the 
numerical examples of this report). 

Transformation of coordinates. — For thin, straight blades 
lying on conic radii (elements), the transformation of co- 
ordinates is given directly by the analytic function 


from which 


and 


^ iv= Ad loge exp 


log e R 
* Ad 




AS 



Figure 9.— Compressor-design characteristics for numerical examples. 


so that, from equations (16a), (16b), and (16c) the coeffi- 
cients in equations (18) and (19) become 

Ui~ 0 


and 


v *-Vt-RAe 

Incompressible solution. — The incompressible solution 
was obtained from equation (18) for the same impeller-tip 
Mach number M T and for the same flow coefficient 0 used 
in the compressible solution but with the density ratio p/p 0 
consant and equal to 1.0. Because for incompressible fluids 
the speed of sound is infinite, M T} 0, and the velocity ratios 
for the incompressible solution are fictitious cpiantities, the 
definitions of which contain a constant, finite speed of sound 
that is equal to c a for the compressible solution. The same 
value of the impeller-tip speed (and of the compressor flow 
rate) therefore results from the same value of M T (and of 0) 
for the compressible and incompressible solutions. 

Design characteristics and operating conditions. — The 
numerical examples have been computed for the following 
design characteristics and operating conditions: 

Design characteristics: 

(1) Constant flow area normal to conic flow surface, 

H=R~ l 

(2) Cone angle a , 180° (or certain other values of a less 

than 180° given by equation (20) for the same value 
of Ad but for different integral values of B) 

(3) Straight thin blades along radii 

(4) Number of blades B } 20 (or other integral values of B 

less than 20 for the same value of Ad but for certain 
different values of a. less than 180° given by equa- 
tion (20)) 

(5) Whir l ratio upstream of blades 0 

(6) Tip Mach number M T , 1.5 

(7) Flow coefficient 0, 0.5 

(8) Ratio of specific heats 7, 1.4 (for compressible solution 

only) 


863977 — 50 3 
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(a) Compressible-flow example. 

Figure 10.— Relative streamlines for flow through centrifugal compressor with straight 
blades. Streamline designation indicates ratio of flow between streamline and positive 
blade surface (right side of passage) to total flow through passage. Angular blade spacing, 
18°; impeller-tip Mach number. 1.5; flow coefficient, 0.5: constant flow area. 

From equation (20) the included angle A6 between blades 
on the conic flow surface is equal to 18°. The results of the 
numerical examples are presented in figures 10 to 12. These 
figures are discussed in the following paragraphs. 

Streamlines, — The streamline configurations (relative to 
the impeller) for the two examples are shown in figure 10. 
The streamlines are designated in such a manner (^/VO that 
the value of a streamline indicates the ratio of the flow that 
lies between the streamline and the positive blade surface to 
the total flow in the passage. For a given density ratio, the 
streamline spacing is indicative of the velocities relative to 
the impeller, with close spacing indicating high velocities 
and wide spacing indicating low velocities. 

In the compressible-flow example (fig. 10(a)), an eddy is 
attached to the positive blade surface. The fluid in this 
eddy rotates (relative to the impeller) in the opposite direc- 
tion to that of the impeller so that the absolute motion of the 
fluid is irrotational. The size of the eddy (for a given im- 
peller) depends on the relative magnitudes of the volume 
flow rate through the compressor and the impeller-tip speed. 
If the flow rate is zero through the rotating impeller, the eddy 
occupies the entire flow passage and as the compressor flow 


Relative 
s treamhnea 



(b) Incomprossible-flow example. (Tip Mach number and flow coefficient for the incom- 
pressible example are fictitious quantities based on the speed of sound Co for compressible 
example.) 

Figure 10— Concluded. Relative streamlines for flow through centrifual compressor with 
straight blades. Streamline designation indicates ratio of flow between streamline and 
positive blade surface (right side of passage) to total flow through passage. Angular blade 
spacing, H°: impeller-tip Mach number, 1.5; flow coefficient, 0.5: constant flow area. 

rate increases (for the same impeller-tip speed) the eddy de- 
creases in size until it finally disappears. The flow rate at 
which the eddy disappears increases as the impeller-tip speed 
increases. The eddy does not exist in the incompressible- 
flow example (fig. 10(b)) because, although the weight flow 
rate is the same for both examples, the volume flow rate is 
higher for the incompressible-flow example as a result of the 
lower fluid density in the region investigated. 

The flow directions in the vaneless diffuser are greatly 
different for the compressible- and incompressible-flow ex- 
amples. (Compare figs. 10(a) and 10(b).) This difference, 
results from the higher volume flowrate for the ineompressible- 
flow example. This higher volume flow rate requires 
higher radial velocities so that for the same tangential 
velocities the flow directions are different in the two examples. 
(From considerations of constant moment of momentum in 
the vaneless diffuser, the tangential velocities should be about 
the same in both examples because the tangential velocities 
are about the same at the impeller tip. See subsequent 
section Lines of constant pressure ratio.) 
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(a) Compressible-flow example. 

Figure 11 .— Lines of constant velocity ratio relative to impeller. 

Lines of constant relative velocity ratio. — Lines of constant 
velocity ratio relative to the impeller are shown for the two 
examples in figure 1 1 . The constant c 0 in the denominator 
of the velocity ratio is the same for both examples and is 
equal to the absolute stagnation speed of sound upstream of 
the impeller for the compressible-flow example. The general 
characteristics of these plots are similar. The velocities 
along the negative blade surface are higher than along the 
positive blade surface except at the tip of the blade where 
the velocities become equal on both blade surfaces (as re- 
quired by the Joukowski condition) . The maximum velocity 
occurs on the negative blade surface at a radius ratio well 
within the impeller and the flow decelerates along the surface 
of the blade from this point to the blade tip. This decelera- 
tion, which becomes rapid near the blade tip, is conducive to 
boundary-layer separation, which lowers the compressor 
efficiency. If the boundary -layer wake in the vaneless dif- 
fuser is neglected, the velocities become essentially uniform 
at a radius ratio of about 1.15. 

In the compressible example (fig. 11(a)), the velocity 
ratios are low at the impeller tip because of the high density 
ratios that result from the high tip speed of the impeller. 
These velocities would be considerably higher if the effective 
flow area were reduced by boundary-layer separation, which 
might be expected in a real compressor. 


CO 


Relative 
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(b) Incompressible-flow example. 

Figure 11 .— Concluded. Lines of constant velocity ratio relative to impeller. 

Lines of constant pressure ratio. — Lines of constant static- 
pressure ratio (local static pressure divided by absolute 
stagnation pressure upstream of the blades) are shown for 
the two examples in figure 12. The general characteristics 
of these plots arc the same. The pressure is higher on the 
positive blade surface than on the negative blade surface 
except at the blade tip where the pressures are equal. This 
difference in pressure accounts for the impeller torque. 

The higher pressure ratios in the compressible-flow example 
than in the incompressible-flow example result from the 
lower relative velocity ratios in the compressible -flow ex- 
ample and from the fact that for the same amount of work 
per pound of fluid the pressure ratio is greater for compressible 
than for incompressible fluids. (That the work per pound 
of fluid is about the same for both examples at correspond- 
ing points is seen from the last term in equation (8). This 
term is the work per pound of fluid and has about the same 
values for both examples because the whirl ratio X is deter- 
mined principally by the tangential motion of the blades, 
which is the same in both examples.) 

Slip factor. — The impeller slip factor is defined as the ratio 
of the average absolute tangential velocity of the fluid at the 
impeller tip to the impeller-tip speed. A method for com- 
puting the slip factor from a relaxation solution is outlined 
in appendix H. The slip factor is 0.899 for the compressible- 
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(a) Compressible-flow example. 

Figure 12 —Lines of constant pressure ratio (local static pressure divided by absolute stagna- 
tion pressure upstream of impeller). 

flow example and 0.892 for the incompressible example. 
It is concluded that the slip factors are essentially equal 
for both examples. 

Compressibility effects. — Figure 10 indicates a large com- 
pressibility effect upon the streamline configuration in high- 
speed, rotating, radial- and mixed-flow blade systems. 
Figures 11 and 12 indicate large compressibility effects 
upon the magnitudes of the velocity ratios and pressure 
ratios, but the distribution of these quantities is similar. 
(For example, the velocities accelerate and decelerate at 
approximately Corresponding positions of the flow field in 
both examples.) It is concluded that, if the fluid in high- 
speed, rotating, radial- and mixed-flow turbomachines is 
compressible, incompressible solutions give poor quanti- 
tative results (exception, the slip factor) and, in some 
respects, poor qualitative results. 

SUMMARY OF RESULTS AND CONCLUSIONS 

A general method of analysis has been developed for two- 
dimensional, steady, compressible flow in stators or rotors of 
radial- and mixed-flow turbomachines with arbitrary blade 
shapes, arbitrary variations in the passage height, and with 


U) 



Angle } deg 


(b) Incompressible-flow example. 

Figure 12. — Concluded. Lines of constant pressure ratio (local static pressure divided by 
absolute stagnation pressure upstream of impeller). 

conic flow surfaces (surfaces of right circular cones generated 
by center line of flow passage in axial-radial plane) . 

The analysis indicates that: (1) The solution obtained for 
a given turbomachine also applies to certain other (equiva- 
lent) turbomachines with a larger or smaller number of like 
passages (same spacing of the blades on the conic flow sur- 
face, same blade-thickness distribution, and so forth) but 
with different cone angles; (2) for the same number of 
similar blades, the blade loading is less for mixed-flow than 
for radial-flow turbomachines; and (3) any solution obtained 
for an outflow turbomachine with shockless (smooth) entry 
is also the solution for an inflow turbomachine with shockless 
entry and with the flow direction and blade rotation (if any) 
reversed. 

Two numerical examples are presented — one for compres- 
sible and the other for incompressible flow in a centrifugal 
compressor with thin, straight blades lying on conic radii 
(elements). The solutions were obtained in a region of the 
compressor, including the impeller tip, that was assumed 
to be unaffected by the inlet configuration of the impeller 
or by the diffuser vanes (if any). Both examples are for 
the same impeller (18° included angle between blades on 
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the conic flow surface) with the same tip speed (equivalent 
to a tip Mach number of 1.5 for the compressible-flow 
example), with the same flow rate (flow coefficient, 0.5), 
and with constant flow area normal to the conic flow surface. 
The following results were obtained: 

(1) In the compressible-flow example, an eddy is attached 
to the positive blade surface. The fluid in this eddy rotates 
in the opposite direction to that of the impeller. This eddy 
does not exist in the incompressible-flow example. 

(2) In both examples, the maximum velocity occurs on 
the negative blade surface at a radius ratio well within the 
impeller and the flow decelerates along the surface of the 
blade from this point to the blade tip. This deceleration, 
which becomes rapid near the blade tip, is conducive to 
boundary-layer separation, which lowers the compressor 
efficiency. 


(3) If the boundary-layer wake in the vaneless portion 
of the diffuser is neglected, the velocities become essentially 
uniform at a radius ratio of about 1.15. 

(4) The slip factor is 0.899 for the compressible-flow 
example and 0.892 for the incompressible example. It is 
concluded that the slip factors are essentially equal for 
both cases. 

(5) If the fluid in high-speed, rotating, radial- and mixed- 
flow turbomachines is compressible, incompressible solutions 
give poor quantitative results (exception, the slip factor) 
and, in some respects, poor qualitative results. 

Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, November 1> 1948 . 



APPENDIX A 


SYMBOLS 


The following symbols are used in this report: 
a T flow area of annulus at tip of blades 

B number of blades (or passages) 

b grid spacing (fig. 3) 

c local speed of sound 

c p specific heat at constant pressure 

e Cartesian coordinate in transformed e/-plane, 

equation (Cl a) 

exp exponential, [exp(x)^e x ] 

F any twice-differentiable function of two variables 

/ Cartesian coordinate in transformed ef- plane, 

equation (Clb) 

g acceleration due to gravity 

H passage-height ratio, h/h T 

h passage height normal to conic flow surface (fig. 1) 

J mechanical equivalent of heat 

M t blade-tip Mach number, equation (5) 

Q relative velocity ratio, q/c 0 

q velocity of fluid relative to blades, yfu 2 +v 2 

R conic-radius ratio (coordinate of conic flow sur- 

face, 220-plane) (fig. 2), r/r T 
R residual 

r conic radius (distance along conic element from 

apex of cone) (fig. 1) 

T static (stream) temperature 

U relative tangential-velocity ratio, ujc 0 (fig. 2) 

u tangential component of q (positive in direction 

of increasing 0) 

V radial-velocity ratio, v/c 0 (fig. 2) 

v radial (along conic element) component of q 

(positive in direction of increasing radius ratio) 
W total flow rate through turbomachine 

w flow rate between streamlines 

a cone angle (fig. 1) 

r dimensionless absolute circulation 

y ratio of specific heats 

A0 angular blade spacing (included angle between 

blade camber lines in 220-plane), equation (20) 
Ai p changes in \p at grid points 

t? Cartesian coordinate in transformed f^-plane 

(corresponds to incompressible stream function 
in nonrotating 220-plane with constant passage 
height, 27=1), equation (13a) 

6 angle (coordinate of conic flow surface, 220-plane) 

(fig- 2) 

X whirl ratio, equation (9) 


£ Cartesian coordinate in transformed ^-plane 

(corresponds to incompressible velocity poten- 
tial in nonrotating 220-plane with constant pas- 
sage height, 72=1), equation (13a) 
p weight density of fluid 

<t> flow coefficient, equation (24) 

\p dimensionless, compressible stream function, 

equations (4a) and (4b) 

co angular velocity of rotor (in direction of increas- 

ing 6) 

Subscripts: 

A,B,C grid points defined in figure 7 

a,b,c,d grid points defined in figure 8 

D downstream boundary (boundary in region of 

uniform flow downstream of blades) (fig. 4) 
i indicates that u iy v if and q i} obtained from deriva- 

tives of £(22,0) and 77 ( 22 , 0 ), correspond to incom- 
pressible velocities 

n negative blade surface (blade surface in direction 

of decreasing 0) (fig. 4(a)) 

o absolute stagnation condition in region of uniform 

flow upstream of blades 

p positive blade surface (blade surface in direction 

of increasing 0) (fig. 4(a)) 

T blade tip (either nose or tail of blade, whichever 

has larger conic radius) 

U upstream boundary (boundary in region of uni- 

form flow upstream of blades) (fig. 4) 

R,6£,r],e, partial derivatives with respect to these coordi- 
and / nates 

RR } 09, second partial derivatives with respect to these 
££, rfrf, coordinates 

£ 77 , ee, 
and f 

1,2, 3, 4 grid points adjacent to point being considered 

(fig- 3) 

Superscripts: 

b' nonuniform grid spacing defined in figure 7 

A change in (at any grid point) resulting from 
At p D required to satisfy Joukowski condition 
A\l/ U change in \f/ (at any gild point) resulting from 

At Pd equal to unity 

adjusted value of \p (at any grid point) after 
Joukowski condition is satisfied 
1 ^* value of (at any point in region of uniform flow 

upstream of blades) if specified values of X^ is 
obtained, equation (F2) 
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APPENDIX B 

TRANSFORMATION OF COORDINATES FROM R0-PLANE TO fr-PLANE 

2 M t H £= (*«+*„) (««*+»,*) +*, (vBB+jf+jp )- t 

h (^B+ff+jp)- 

(.'hVt—'l'vU’t)[(\oge H)iVi~ (loge H\U t ] — 

[+• ( ,og ‘ £),+*' ( log -£),] w+vf> 


If the transformation of coordinates from the i?0-plane to 
the ^-plane is represented by the analytic function 


exp (id)] 


(13) 


where the coordinates £ and rj in the £^-plane correspond to 
velocity potential lines (£= constant) and streamlines (tj~ 
constant) in the 7?0-plane for incompressible flow past the 
blades, which for purposes of the transformation are con- 
sidered to be stationary (<*>=0) and to have a constant height 
(H—l), then 


Ui=~r} R = 


& 

R 


V 9 y 

Vi—R — tR 


(16a) 

(16b) 


Also, if F is any twice-differentiable function of R and Q 


F r = F£ r + F v 7] r 

F rr = F^r 2 ~r 2 F$ v £ r t} R + F vv rj R 2 + F £ RR + F v rj RR 
F 6 = F£q + F v 7]0 

Fee —Ff.&e 2J r2F%r} Vo + F vv r} 0 2 + Ft £ e9 + F^ee 


(Bl) 


From equations (16a), (16b), and (Bl), equation (7) 
becomes 


But, 

and, because equation (13) is analytic, 


and 


£rr + %+R2—° 


^/e + ^ + ^2^0 


so that equation (7) finally becomes 

^*r *> ( l0 *v,)r * ( log - Or 

tPi-pV'j [(loge (]oge h\u J 

Equation (12) in like manner becomes 

«£-hOw+'M“ 


(16c) 


(14) 


(15) 
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APPENDIX C 

${R 9 9) AND v(R,0) for arbitrary blade shapes 


The coordinates £ and 77 of the transformed £77 -plane are 
functions of R and 0 and correspond to the velocity potential 
£(I2,0) and the stream function 77(JR,0) for incompressible flow 
past blades of arbitrary shape in the physical l?0-plane, 
which blades are considered, for purposes of the transforma- 
tion, to be stationary (w=0) and to have a constant height 
(. H=l ). In order to determine £{R,Q) and ??(1?,0) ; it is con- 
venient first to transform the blades from the I? 0-plane to the 
e/-plane (ef Cartesian coordinates). This transformation is 
given by 

e+if= log e [R exp (iff)] 

from which 

e—lo g e R (Cl a) 

/=0 (Clb) 

Equations (Cla) and (Clb) relate points in the ^/-plane to 
points in the l?0-plane and determine a new blade shape in 
the ef - plane (fig. 13) that corresponds to the original (arbi- 
trary) blade shape in the 7?0-plane (fig. 4(a)). In effect the 
radial cascade in the R 0-plane is transformed into an axial 
cascade in the ^/-plane. 



The stream function rj(ef) in the <?/-plane is determined by 
the relaxation solution of Laplace’s equation 

Vee+V/f—Q (C2) 

for the specified boundary conditions; and the velocity 
potential £(e,j0 is obtained from 7j(e,f) by methods discussed 
in reference 7 (ch. IV). Finally, £(2?,0) and rj(R y 6) are de- 
termined from £(e,/) and and from equations (Cla) and 

(Clb). 

In order to solve equation (C2), it is first necessary to 
determine the boundary conditions. Because the distribu- 
tion of the variations in v along lines of constant e in the 
ef- plane is cyclic with a period equal to the blade spacing, 
the solution of equation (C2) need be obtained only in a 
region (fig. 13) bounded by the surfaces of two adjacent 
blades, by lines of constant e that correspond in the i?0-plane to 
the upstream and downstream boundaries of the compressible- 
flow field (fig. 4 (a)), and by quasi boundaries extending along 
lines of constant j from the ends of the blade surfaces to the 
upstream and downstream boundaries. (Note that these 
quasi boundaries do not generally, and need not. correspond 
to the quasi boundaries selected for the compressible-flow 
field (lines of constant 77). See fig. 13.) 

After the location of the boundaries has been determined, 
a grid of equally spaced points is placed inside the boundaries 
(fig. 13) and this grid is extended to points on the boundaries. 
The grid points on the blade boundaries generally are un- 
equally spaced from the interior grid points because of the 
arbitrary shape of the blades. Values of 77 at points along 
the quasi boundaries extending from the negative blade 
surface are directly related to values of 77 at points along the 
quasi boundaries extending from the positive blade surface 
(see related discussion in section ^ along quasi boundaries 
extending from positive blade surface) and, therefore, only 
the values of 77 along the quasi boundaries extending from the 
positive blade surface need be recorded (fig. 13) and relaxed. 

The values of 77 at the grid points on the boundaries are 
next determined. The boundary values of 77 at grid points 
along the blade surface are arbitrarily set equal to 0 along 
the positive blade surface (fig. 13) and equal to 1.0 along the 
negative blade surface. The values of 77 at grid points along 
the quasi boundaries extending from the positive blade surface 
are estimated in such a manner that the front and rear 
stagnation points occur at the intersection of the mean 
camber line with the surface of the blade at the nose and the 
tail. The direction of the streamlines is then approximately 
equal to the direction of the blade camber line at the ends 
of the blade. This direction is defined in the e/-plane by 

ten p=f g 


Figure 13.— Relaxation grid in e/-plane used to obtain y(ej). 


(C3) 
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But along a streamline, dij equals zero so that 


dt]=0=i]jle+y J df 

if = _ve 

de 7jf 


(04) 


Also, if uniform flow is assumed upstream and downstream 
of the blades, then from equation (Clb) and the specified 
boundary conditions 


Vf Af A 6 

Therefore, from equations (C3), (C4), and (C 5 ), 


(05) 




tan /? 
Ad 


( 06 ) 


This variation in 77 with e determines estimated values of 77 
at grid points along the quasi boundaries extending from the 
positive blade surface and in particular this variation deter- 
mines the values of 77 at the intersections of these quasi 
boundaries with the upstream and downstream boundaries 
(fig. 13). Along these upstream and downstream boundaries, 
the values of 77 increase uniformly (steady-flow condition) 
in the positive direction of / at a unit rate per blade spacing. 

After the boundary values of 77 have been determined, 
values of 77 are estimated at the interior grid points. These 
estimated values are generally in error and must be corrected 
by relaxation methods in which equation (C 2 ) is used in 
finite-difference form to compute and to relax the residuals. 
In addition, the values of 77 along the quasi boundaries extend- 
ing from the positive blade surface were estimated values 
and must therefore be relaxed. After the solution for the 
distribution of 77 has been obtained, the condition that the 
stagnation points occur at the nose and the tail of the blades is 


checked (by methods similar to those in appendix E) and, 
if not satisfied, the values of 77 at the grid points along the 
upstream and downstream boundaries are adjusted (by 
methods similar to those outlined in appendix E). 

The function 77 (e,f) is now known and £(ej) can be deter- 
mined from 77 (e,f) by methods given in reference 7 (ch. IV). 

The functions 77 (R,6) and £(i?,0) are obtained directly 
from y(ej) and £(e,f) and from equations (Cl a) and (Clb). 
Also, from equations (16) and (Cl), the coefficients u t and 
Vi in equations (18) and (19) become 

— rjR= — ^ (C7a) 


Vi= 


ye_yf 

R~R 


(07b) 


The relaxation solution of equations (11), (14), and (15) 
in the transformed ^ 77 -plane requires less time than the solu- 
tion of equations (7), (11), and (12) in the physical A?0-plane, 
because blades with arbitrary shape in the physical plane 
become straight and parallel in the transformed plane, 
which results in simpler finite-difference forms for equa- 
tions (14) and (15). The transformation of coordinates to the 
^ 77 -plane is time-consuming, however, so that, if a solution 
for only one set of operating conditions is desired, it would 
probably be faster to solve equations (7), (11), and (12) in 
the e/-plane where, although the finite-difference equations 
must contain coefficients to account for the unequal grid 
spacing along the irregular boundaries, the transformation 
of coordinates is given directly by equations (Cl a) and (Clb) . 
If, however, solutions for a number of different operating 
conditions for the same blade configuration are desired, then 
the transformation of coordinates outlined in this appendix 
is desirable, because the same transformation applies to all 
sets of operating conditions for the same blade configuration. 


APPENDIX D 


ESTIMATED VALUES OF ^ AT GRID POINTS ALONG QUASI BOUNDARIES EXTENDING FROM POSITIVE BLADE SURFACE 

IN £t7-PLANE 


Estimated values of the stream function $ at grid points 
along the quasi boundaries extending from the positive blade 
surface in the ^-plane can be obtained by assuming, as a 
first approximation, that the flow conditions upstream and 
downstream of the blades in the i20-plane are uniform, that 
is, the flow conditions are a function of R only. From the 
conservation of absolute moment of momentum (whirl) 
upstream and downstream of the blades 


where, from equation (3a), 

Q2=U 2 +V 2 (D5) 

But, from continuity considerations assuming uniform flow, 



X = R ( RM t + U) = cons tan t (Dl) 

so that equation (4b) becomes 

* r =^h(rM t -j^ (D2) 

In addition, because the flow is considered uniform 

ype = constant = (D3) 

The variation in \p along the quasi boundaries in the f??-plane 
(fig. 6) is then given by 

which, from equations (D2) and (D3) becomes 

H (^RM t —^J 0* CD4) 

where R * and 9% are obtained from equations (13a) or ap- 
pendix C. Equation (D4) gives the estimated variation in 
\p along the quasi boundaries extending from the positive 
blade surface in the £17- plane. 

In order to integrate equation (D4), it is necessary to 
know the variation in density with £. The density ratio is 
given by equation (11) 

i 

p P o ={ 1 + 7 ~ 1 [( RM T )' 2 — Q 2 — 2 M r Xr /]| T_1 ( 11 ) 


so that, from equations (Dl), (D5), and (D6), 



and equation (11) becomes 



(D7) 


Because R is a known function of £ and tj (equation (13a) or 
appendix C), the system of equations (D4) and (D7) can be 
solved by numerical methods to obtain the value of \p at 
grid points along the quasi boundaries extending from the 
positive blade surface in the £??-plane (fig. 6). 

The values of \p depend on the value of X. Upstream of 
the blades X has the specified value \ v . Downstream of the 
blades X has the value \ D , which for a given blade shape and 
operating condition, is determined by the Joukowski condi- 
tion. As a result of the Joukowski condition, the average 
flow direction at the exit from the blades is approximately 
equal to the blade-exit angle (determined by the mean 
camber line) . An average value of U, required in equation (9) 
to compute an estimated value of \ D , can therefore be 
obtained from this angle (adjusted as experience indicates) 
and from the average value of V given by continuity con- 
siderations (equation (D6)). 
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APPENDIX E 

METHOD OF ADJUSTING VALUES OF $ ALONG DOWNSTREAM BOUNDARY TO SATISFY JOUKOWSKI CONDITION 


The Joukowski condition requires the rear stagnation point 
to occur at the blade tail, or, in case of infinitely thin blades 
or blades with cusped tails, the flow must be tangent to the 
blade surfaces at the tail. (If the blade tail is somewhat 
rounded, the stagnation point is considered to occur at the 
intersection of the mean camber line with the tail surface 
of the blade.) In the f??-plane (fig. 6), this rear stagnation 
point occurs at the tail of the thin, straight blade. (See 
appendix C.) This condition generally is not satisfied by 
the initial relaxation solution for \p in the ij^-plane, because 
for this solution the values of \p along the downstream bound- 
ary (i p D ) were obtained from the estimated variation in \p 
along the quasi boundaries (appendix D) . In order to satisfy 
the Joukowski condition, the values of \p at the grid points 
along the downstream boundary (\p D ) must all be changed 
the same required amount (A \p D ). This change in \p D 
(denoted by A\p D ) results in changes in \p (denoted by Axp) 
at each of the interior grid points and at the grid points 
along the quasi boundaries. The manner in which the values 
of \p are changed by the change in \p D must satisfy the differ- 
ence equation (18). Therefore, 


(yp\ + Ai/'J + (\p 2 + A \p 2 ) + (^3 + A^ 3 ) + (1 p.i + Ai p 4 ) — 4 (\p + A xp) — 


W'i +_A ^j) — i)p‘i + A^ 3 ) 
4 

GA2 + - A1A2) — (\ A4 + A 1^4) 
4 


( 


Pi _ 
Po 


r *Y 

' Po) 

(log. *-1og. *)- 

\ Po Po/ 


The solution of equation (Elb) determines A \p at every 
grid point for a specified value of A \p D . Because of the 
linearity of equation (Elb), the solution for any specified 
value of A\p D is equal to the solution for A\p D = 1.0 multiplied 
by the specified value of A\p n . That is, A\p (at any grid 
point) resulting from a specified value of Ai p D is equal to Ai p 
(at the grid point) resulting for A\p D =l.Q multiplied by the 
specified value of A >p D . 

The procedure for the solution of equation (Elb) is the 
same as for equation (18). The boundary values of \p along 
the blade surfaces and along the upstream boundary (fig. 6) 
are not changed so that A \p must equal zero along these 
boundaries. The value of A \p along the downstream bound- 
ary (fig. 6) is set equal to unity. 

The magnitude of A \p D required to satisfy the Joukowski 
condition can now be determined as follows: If the rear 
stagnation point occurs at the blade tail (Joukowski condi- 
tion), then the extrapolated value for \p at the grid point on 
the blade tail of the positive blade surface obtained from 
the values of \p at succeeding points along the quasi boundary 
starting at the blade tail must equal zero. The extrapolated 
value of \p at the blade tail using a third-degree polynomial 
and the first four points along the quasi boundary is given by 

'Pluil = 4'Pa — 6'Pb+4\p c —'ipd 

where the subscripts a, b , c, and d refer to the grid points 
along the quasi boundary in figure 8. If \p fa u equals zero, 


4 1 2 {[(*.+ A*,) - (*, + A *,)]»,- [ (*,+ A* 2 )-(* 4 + A* 4 ) }u t } X 


{ (log. III— log. II 3 )>'i— (log,. Hi — log, III)Ut} — 
2 M T rib 2 P = R 

n / o 


(El) 


where the change in density ratio resulting from A \p is con- 
sidered negligible. Subtracting equation (18) with R equal 
to zero (which condition has been satisfied by the initial 
relaxation) from equation (El) results in 


0=4*.‘-6* ft < +4* 6 *-*,‘ (E2) 

where \p l signifies values of xp after the Joukowski condition 
is satisfied. But, 

yp i =yp J r Axp 1 (E3) 

where \p is the stream function obtained by the initial 
relaxation and Ai p l is the change in \p that results when the 
Joukowski condition is satisfied. Also, from the first part 
of this appendix, 

A^=A^ d A^* (E4) 


Af. + Afc + Afc + Afc-^A*- 4 * 1 A * 3 (log. Pl -log. p3 )- 

Ati-Mu P2_ loge P ^) 4^- 2 [(A^j Atf>3)Vi (A^2 -A tju,] X 

[ (log. Hi — log. II 3 ) Vi— (log. Ih — log. «,] = R (E 1 a) 

Each of the last three terms on the left side of equal t.ion (El a) 
consists of the product of two quantities that approach zero 
as the grid spacing b approaches zero. For the small grid 
spacing used in relaxation solutions, these terms are there- 
fore of secondary importance and may be neglected so that 

Axp 1 -\-A\p 2 + A^ 3 + A^ 4 — 4A\p~ R (Elb) 


where A\p li is the change in \p (at any grid point) resulting 
from a unit change in \p D (Axp D =l.O) and Axp D is the change in 
\p D required to satisfy the Joukowski condition. Therefore, 
from equations (E2), (E3), and (E4)) 


A , "4^ g +6^“4^+^ 

4 A^ a * * — 6 A^ & if + 4 A^J * — A\p d u 


(E5) 


Equation (E5) determines the change in \p D required to 
satisfy the Joukowski condition. The changes in \p at all 
other grid points are obtained by multiplying A \p H at each 
grid point by A \p D . Because the solution for A \p H is approxi- 
mate, the resulting values of 1 p l must usually be relaxed to 
eliminate small residuals computed by equation (18). 
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APPENDIX F 

METHOD OF ADJUSTING VALUES OF * ALONG UPSTREAM BOUNDARY TO OBTAIN SPECIFIED WHIRL RATIO X^ 


In general, the specified whirl ratio upstream of the blades 
Xu is not obtained by the initial relaxation solution, because 
for this solution the values of xp along the upstream boundary 
(xp v ) were obtained from the estimated variation in xp along 
the quasi boundaries (appendix D). In order to obtain the 
specified value of Xu, the values of xp at the grid points along 
the upstream boundary (xp v ) must all be changed the same 
required amount (A^). This change in xp v (denoted by 
A xpu) results in changes in xp (denoted by At p) at each of the 
interior grid points and at the grid points along the quasi 
boundaries. The effect of A xp v on the values of Ai p is deter- 
mined in the same manner as the effect of Axp D on the values 
of At p (appendix E) . 

The magnitude of A xp v required to obtain the specified 
value of X v can now be determined as follows: Near the 
upstream boundary in the region where flow conditions are 
essentially uniform, the whirl ratio X is constant and equal 
to Xc/. In this region equation (9) gives 

\u=R U') (9) 

where U is related to the variation in xp by equation (4b) 


But, 





where for uniform flow conditions 




(4b) 

(Fl) 


so that, from equations (16a) and (16b), equation (Fl) 
becomes 

'pR='PiV i —\p rl Ui (Fla) 
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Therefore, from equations (Fla), (4b), and (9) 

(F2) 


where xp e * is the value of xp$ at any point in the region of 
uniform flow upstream of the blades if the specified whirl 
ratio \u is obtained and where p/p 0 is considered to be given 
by the initial relaxation solution. (In general, equation (F2) 
is evaluated at the upstream boundary where, assuming- 
conditions are uniform, xp% is constant.) 

If xp$ is obtained by the initial relaxation solution at the 
point being considered in the region of uniform flow, the A \p v 
required to obtain the value of x p%* (equation (F2)) corre- 
sponding to the specified value of Xu is given by 

(F3) 


where (Ai p)z is the variation in Ai p with £, at the point being 
considered, for a unit change in x p v . Equation (F3) deter- 
mines the change in x p v required to obtain the specified value 
of X v . 

The resulting changes in xp at the interior grid points and 
at the grid points along the quasi boundaries are determined 
from Axpu in the same manner as the changes in xp were deter- 
mined from Ai p D in appendix E. It should be noted that the 
correction for Xu will affect the Joukowski condition, and vice 
versa. For high-solidity blades, however, the effect of Axp D 
on X v and the effect of Ai p v on the Joukowski condition is 
generally small and can be neglected. 
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APPENDIX G 


SIMPLIFIED ANALYSIS FOR ROTORS WITH STRAIGHT BLADES ALONG CONIC RADII 


The relaxation methods used in this report are lengthy. 
It would therefore be advantageous to have a quicker, al- 
though less accurate, means of estimating the flow conditions. 
In this appendix, a simplified analysis is developed for rotors 
with straight blades along conic radii. 

Velocity-ratio distribution. — This simplified analysis is 
based on the assumption that for rotors with straight blades 
along conic radii the tangential component of the velocity 
ratio relative to the blades is zero within the rotor. Equa- 
tion (6) therefore reduces to 

V a =2 BM t (Gl) 

which, when integrated, becomes 

V=V P +2 RM t 6 (G2) 


where the subscript p refers to the positive blade surface at 
which surface the angle 6 is considered zero. Equation (G2) 
gives the distribution of the radial component of the velocity 
ratio across the passage at constant values of R. The con- 
stant of integration V v is determined at each value of R 
from considerations given in the next paragraph. 

Stream-function distribution. — From continuity 


or 


dw= pvhrdd 


dw 

PoCoh-^T t 


L VHRdd 

Po 


(G3) 


The density ratio is given by equation (11) with Q equal to 
V (because U is assumed equal to zero) and V is given by 
equation (G2) so that equation (G3) becomes 

d*=HR jl+^p [{RM T y-(V p + 

1 

2RM t 6) 2 2A/r\t/]j ( V p -\- 2RM T d)d6 (G4) 


where the left side of equation (G4) was obtained from 
equation (21). Equation (G4) is integrated from the posi- 
tive blade surface where and d are considered equal to 
zero so that 

' p= mr T \_\ 1+T ^ 1 mM T y-vy- 2 M T \ v ]^~ l - 
{(RM T y- (V p + 2 RM T ey- 2 M T \ v ] p 1 ] 

(G5) 

Equation (G5) gives the distribution of across the passage 


at constant values of R. The velocity ratio V p varies with 
R and is obtained from equation (G5) for the condition 


when 




6 — 6 n 


If the fluid is incompressible, the distribution of \p becomes 

t=HR(V p d+RM T d 2 ) 

Numerical example. — The simplified analysis has been 
applied to the compressible-flow example in this report and 
the results are compared with those of the relaxation solu- 
tion. The velocity ratio V p along the positive blade surface 
has been computed from equation (G5) and the results are 
compared in figure 14 with the relaxation solution. The 
negative values of V p occur where the eddy (fig. 10(a)) is 
attached to the blade. The agreement between the re- 
laxation solution and the simplified solution is satisfactory 
up to a radius ratio of about 0.80. For radius ratios greater 
than 0.80, the agreement is unsatisfactory because the 
assumption that U and its derivative are negligible is no 
longer valid. 



Figure 14. — Comparison of velocity ratio along positive blade surface for simplified and 
relaxation solutions. 
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Ratio of distance across passage 
from positive blade surface 

Figure 15 —Comparison of velocity distribution across passage for simplified and relaxation 

solutions. 

The velocity-ratio distributions across the passage at radius 
ratios of 0.760 and 0.855 have been computed from equation. 
(G2) using values of V p obtained from equation (G5) and 
the results are compared in figure 15 with the relaxation 
solution. At the 0.760 radius ratio, the velocity distribution 
is nearly the same for both solutions, but at the 0.855 radius 
ratio the simplified solution has begun to deviate from the 
more rigorous relaxation solution. 



Figure 16.— Comparison of stream-function distribution across passage for simplified and 

relaxation solutions. 


The stream-function distributions across the passage at 
radius ratios of 0.760 and 0.855 have been computed from 
equation (G5) and the results are compared in figure 16 
witli the relaxation solution. At a radius ratio of 0.760 the 
stream-function distribution is nearly the same for both 
solutions, but at 0.855 the simplified solution has begun to 
deviate appreciably from the relaxation solution. 





APPENDIX H 

PROCEDURE FOR COMPUTING IMPELLER SLIP FACTOR 


The impeller slip factor for centrifugal compressors is 
defined as the ratio of the average absolute tangential veloc- 
ity of the air at the impeller tip to the tip speed of the 
impeller 


Slip factor = 


sin f+“)„ ,,u, 

I T H A 


c or T sin 2 


+w T (H1) 


The average value of the tangential-velocity ratio relative ' 
to the impeller at the impeller tip is obtained from 


u Q 


or 


U av = 


( n Upvh T r T dd 

j== 

j p'vh'pV'p do 
PoC^t^t C^ n jjy P 

J 0 P Po 


w 

B 


which, from equation (22), becomes 

U av =] f ""uv-^de 

WnjK Po 


(H2) 


Equation (H2) gives a weighted average value of U. 
This weighted average is also equal to the unweighted average 


1 CK 
3 n Op Jo p 


Udo 


This fact can be shown from considerations of the conserva- 
tion of moment of momentum in the vaneless diffuser, which 


is based upon the weighted average of U and from considera- 
tions of constant absolute circulation in the diffuser, which is 
based upon the unweighted average of U. 

Combining equations (Hi) and (H2) results in the follow- 
ing expression for the slip factor : 

Slip factor = 1 P“ UV ~ de (S3) 

The value of the integrand is obtained from the relaxation 
solution. 
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5. NUMERICAL RELATIONS 


- 1 hp=76.04 fcg-m/s=550 ft-lb/sec ,, 
1 metr£ horsepower— 0.9863 Ep 
1 mph= 0.4470 mp 3 V- V A- j 
A 1 mps= 2.2369 mph ■ - . X 


1 lb=0.4536 kg 
1 kg=2.2046 lb ; , 

I mi= 1,609.35 m= 5,280 ft: 
1 m= 3.2808 ft - .... 



